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We analyse forward-jet production at HERA in the framework of the Golec-Biernat and Wiisthoff 
saturation models. We obtain a good description of the forward-jet cross-sections measured by 
the HI and ZEUS collaborations in the two-hard-scale region (kr ~ Q > A-qcd) with two dif- 
ferent parametrizations with either significant or weak saturation effects. The weak saturation 
parametrization gives a scale compatible with the one found for the proton structure function F2. 
We argue that Mueller-Navelet jets at the Tevatron and the LHC could help distinguishing between 
both options. 
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The saturation regime describes the high-density phase of partons in pcrturbative QCD. It may occur for instance 
when the Balitsky-Fadin-Kuraev-Lipatov (BFKL) QCD evolution equation [| goes beyond some energy related to the 
unitarity limit IE El IE IE ■ O n a phenomenological ground, a well-known saturation model [j| by Golec-Biernat 
and Wiisthoff (GBW) gives a parametrization of the proton structure functions already in the HERA energy range. 
It provides a simple and elegant formulation of the transition to saturation. However, there does not yet exist a 
clear confirmation of saturation since the same data can well be explained within the conventional perturbative QCD 
framework 0. 

. In fact, the study of the proton structure functions is a one-hard-scale analysis since their QCD properties are 
dominated by the evolution from a soft scale (the proton scale) to the hard scale of deep inelastic scattering. In 
t-H , order to favor the evolution at fixed transverse scale, which is expected to lead more directly to saturation, it seems 
t-H ■ interesting to focus on two-hard-scale processes such as forward-jet production at HERA, in the region where the 
I transverse momentum of the jet kr is of the order of the virtuality of the photon Q 2 and where the rapidity interval 
for soft radiation is large. This kinematical configuration which was already proposed 10] for testing the BFKL 
evolution is thus also a good testing ground for saturation. The goal of our paper is then to formulate and study the 
. extension of the GBW model to forward jets at HERA. 

Among the physics questions that we want to raise, we analyse whether saturation effects can be sizeable in forward 
jets at HERA and whether they are compatible or not with the GBW parametrization of F^- More generally, we 
would like to compare potential saturation effects in one-scale and two-scale processes. Indeed, in the two-scale 
Q_i! process initiated by 7*7* scattering, it has been suggested 0] that the saturation scale could be different and in fact 
quite larger than the one of deep inelastic scattering. However, another study of 7*7* scattering sticks to the same 
saturation scale as for scattering ^^|. We will discuss this point with forward-jet data which have better statistics 
^ . and wider kinematical range than 7*7* scattering. If the saturation scale is universal as assumed in |12|. then one 
should not expect to see stronger saturation effects in forward-jet data than for F%\ however if the saturation scale is 
higher for processes initiated by two hard scales as proposed in jl l| . saturation effects could be more important. 

For HERA, the existence of two different saturation scales for one or two hard-scale processes would be an interesting 
feature. This study is also of interest in the prospect of Mueller-Navelet jets ^E] a t the Tevatron and the LHC where 
saturation with two hard scales could play a bigger role. 

The plan of the paper is the following. In Section 2, we formulate the extension of the GBW model for forward 
jets. In Section 3, we present the fitting method and give the results of the fits together with several tests of the 
two saturation solutions that we obtain. Section 4 is devoted to a discussion of our results and to predictions for 
Mueller-Navelet jets at hadron colliders that could discriminate between both solutions. 
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FIG. 1: Forward-jet production in partonic and dipole representations. Fig la: Forward jets at HERA. Fig lb: Forward jet in 
7* —onium scattering in the partonic representation (at 1st order, i.e. with only one gluon exchanged). Fig lc: Forward jet 
in 7* — onium scattering in the dipole representation. Y : rapidity gap between the two hard probes, kr > Q : jet tranverse 
momentum lower bound. Q : virtuality of the photon. The gluon-dipole couplings are sketched by black points. 



II. FORMULATION 



The original GBW model provides a simple formulation of the dipole-proton cross-section in terms of a saturation 
scale Rq(Y) ~ e~ 2-( y ~ Y °) where y = logl/ai is the total rapidity, A is the intercept and Yq sets the absolute value. 
Varying Q, QRq 3> 1 corresponds to the dilute limit while when QRq < 1, the dipole-proton cross-section saturates 
to a finite limit a®. In two-hard-scale problems, one has now to deal with dipole-dipole cross-sections which require 
an extension ^1 °f the GBW parametrization. In the case of forward jets, one has to combine the dipole-dipole 
cross-section with an appropriate definition of the coupling of this cross-section to the forward jet. This has been 
proposed in 0] . We shall use this coupling in the analysis of forward jets at HERA. 

The QCD cross-section for forward-jet production in a lepton-proton collision reads 

d [A) cr _ _a_ 2 / / da T da L \ da T y K , 

dxdQHxjdk\ ~ nxQ 2 /e ^ xj ' Kt > \{dk 2 + dk 2 ) 1 V) + dk 2 T o ' ' 



where x and y are the usual kinematic variables of deep inelastic scattering, Q 2 is the virtuality of the photon with 
longitudinal (L) and transverse (T) polarization, and xj is the jet longitudinal momentum fraction with respect to the 
proton. duT.hldk\ are the photon-parton hard differential cross-sections for the production of a forward (gluon) jet 
with transverse momentum fc<r3> h-QCD- The effective parton distribution function / e // has the following expression 

f e f f (xj,kr) = g(xj,kr) + ^(q(xj,kr) + qixj,^)) , (2) 

where g (resp. q, q) are the gluon (resp. quark, antiquark) structure functions in the incident proton. k\ is chosen as 
the QCD factorization scale. 

Properly speaking, the forward-jet cross-section (1) is the leading log(l/a;) expression when one uses the BFKL 
formulation of d<jT,hl dk\. Let us show how one can extend this formalism when saturation corrections are present, 
e.q. in the framework of the GBW model. The differential hard cross-section reads 



do~T,L da 



dk z T dQ 2 



(3) 



where ctt,l are the cross-sections for the production of a forward jet with transverse momentum larger than Q, 



expressed in the dipole framework by [l 



a T , L (Q 2 ,Q 2 ) = ^ j >r d 2 f ^ L (r,Q 2 ) ^(r, Q^) |l - exp ^-^^^ U . (4) 
<7o times the term in brackets in formula @ is the GBW dipole-dipole cross-section. 

Rl(y) = -k e- KY - Yo) (5) 

is the saturation radius redefined for the forward jet case, where Y = \ogxj /x is the rapidity interval available for 
the forward-jet cross-section (see Fig 1), A is the saturation intercept and Y"o is a parameter defining the rapidity for 
which Rq(Yq) = 1/Qo = 1 GeV^ 1 . <tq is the value at which the dipole-dipole cross-section saturates. The effective 
radius r^tAr, f) is defined through the elementary dipole-dipole cross-section given by the two-gluon exchange 0] 

<r dd (r,r) = J d 2 kf°(k 2 ,r)f(k 2 ,f) = 27ra 2 r 2 (l + In ^) = 2nc£ti ff (r,f) , (6) 

where j~. (resp. ri) is min (r, f) (resp. max (r, f)), a s is QCD the coupling constant and f°(k 2 , r) = 2a s (l~ Jo(kr)) / k 2 
is the dipole-gluon coupling. The functions <j>j, L (r, Q 2 ) and (j> J {f, Q 2 ) express the couplings to the virtual photon and 
to the jet. (ftp L (r, Q 2 ) is the known squared QED wavefunction [Tfij | of a photon in terms of qq. 

In equation (0J , the definition of (f> J (f , Q 2 ) associated with a forward jet with transverse momentum kx>Q requires 
more care. The gluon jet is colored while the associated dipole f is not, see Fig.lc. 4> J is obtained by taking advantage 
of the equivalence between the partonic and dipole formulations of forward-jet production when the jet is emitted off 
an onium (see Fig lb-lc, where the size of the incident onium is denoted ro). Let us sketch the derivation of <f> J (f, Q 2 )- 
Assuming the condition 1 GeV^ 1 3>ro 3> 1/Q, the onium is small enough to allow for a perturbative QCD calculation 
but large enough with respect to the inverse transverse momentum of the forward jet. At lowest pertubative order, 
the coupling of the system onium-jet to the gluon (interacting with the target photon) reads (14L \YA Il8j : 

a s log— / ^0(4-Q 2 )/°(|£+fc T |Vo)«(2a s log — logQr l f d 2 f Jx(Qf) f(k 2 ,f) . (7) 
xj J nk\ I xj J J 2nr 

The function f°(k 2 ,f) factorizes as part of the dipole-dipole cross-section ©. The factor in brackets in Q is the 
first order contribution of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) gluon ladder 0, at the Double 
Leading Log (DLL) approximation. QCD factorization thus implies that it is included in the structure function of the 
incident particule. In the case of a proton it is absorbed in f e ff, see formula ©, by a redefinition of the factorization 
scale. Therefore, the remaining term gives 

^(f t Q») = A. ,h(Qf) . (8) 
2irr 



Note that (f> J is an effective distribution that comes out from the calculation rather than a probability distribution as 
Inserting in formula the known Mellin transforms <P(t) = J d 2 r (r 2 Q 2 ) 1_T 0(r, Q 2 ) : 



for <\? T L 



y M -^ r(2-T) ~ 2N c a^ i r 2 (i + 7 )r 2 (i- 7) r 2 (2- 7) /(i + 7 )(2- 7 )\ 

T(r) ' <Pt,lV/)- n Z. e ? 4 7 7 r(2-2 7 )r(2 + 2 7 )(3- 2 7 ) ^ 2 7 (1 - 7 ) ) W 

and after straightforward transformations we can express our results in a double Mellin-transform representation: 

< i?e(r), Re(i), i?e( 7 - r) < 1 , 



where 



a(r, 7 )=a o y du 2 J d* 2 ***-**-*-* (l - e^tt = <r ^ _^ du u^ 1 [r 2 eff (l, u)] 1+T ~ 7 



(11) 



Using JBJ), one finally gets 

a(r, 7 ) = ao — ^ ^ {* 1, 3 + r - 7 , 2r + * 1, 3 + r - 7 , 2 - 2 7 } , (12) 

1 + t - 7 

where the confluent hypergeometric function of Tricomi a, b) can be expressed |Tflj in terms of incomplete Gamma 
functions. 

Note that, expanding the exponential function in the GBW dipole-dipole cross-section in (0J, it is possible to single 
out order by order the contribution of r 2 ™* in formula (|10|1 : 

/_i\n-i 

cr(«)( r 7 ) = 2o- - i- 2i7T(S(7 - r + n - 1) {*(1, 3 + r - 7, 2r) + 3 + r - 7. 2 - 2 7 )} . (13) 

n! 

There is an interesting physical interpretation of this expansion, corresponds to the dilute limit of the GBW 
formula. Taking also into account can be interpreted as saturation a la Gribov-Levin-Ryskin [21, i.e. ~ r 2 ^ — 
r^ff /8Rq(Y). Note also that formulae (|10ll 1|) can also be used with the other models for proposed in the 
literature H3| . 



III. FITTING FORWARD JETS 



For the fitting procedure, we will use a method allowing for a direct comparison of the data with theoretical 
predictions. This method has already been applied [2(| for a BFKL parametrization of the same da/dx data at 
HERA. We shall extend this method to the GBW parametrization (cf. Eq.JTJU)- 

The published data depend on kinematical cuts (see |2lj,|22j) which are modeled by bin-per-bin correction factors 
that multiply the theoretical cross-sections. The details of the method are as follows: i) for each x-bin, one determines 
[23| the average values of Q 2 , xj, kx from a reliable Monte-Carlo simulation of the cross-sections using the Ariadne 
program |24|; ii) one chooses a set of integration variables for d^a (see in such a way to match closely the 
experimental cuts and minimize the variation of the cross-sections over the bin size, the convenient choice of bins for 
forward jets (2£j is 

Ml) <14) 

(note the choice of the variable k\ / Q 2 which is well-suited for the study of two-scale processes) ; iii) one fixes the 
correction factors due to the experimental cuts for each x-bin by a random simulation of the kinematic constraints 
with no dynamical input. The list of correction factors for the HI and ZEUS sets of data [2HI22I and more details on 
the method are given in [2£j ■ 

Using these correction factors, we perform a fit of formulae 0) to the HI and ZEUS data using the GBW 
parametrization (|10I12I) . The free parameters arc the saturation scale parameters A and Y"o (sec (j5Jl) and the nor- 
malizations which we keep independent for HI and ZEUS. Note that they are related to the dimensionless factors 
o~oQ 2 - We also make the slight modification Co ^ a'glk 2 -,) cto with as running at one loop with 4 active flavours and 
Aqcd = 220 MeV. Indeed, it is known that leading-order BFKL fits are better when the coupling constant in the 
overall factor is running. It turns out that with our GBW parametrization, the effect is much smaller than in the 
BFKL case as will be discussed later on. The obtained values of the parameters and the \ 2 of the fits are given in Table 
I and the resulting cross-sections displayed 1 in Fig. 2. We obtain two different solutions which can be characterised 
by different strength of saturation effects. The x 2 per do/ is quite good (we use only statistical errors to perform the 
fit) in both cases. The two solutions show similar \ 2 values and resulting cross-sections and we chose to plot only the 
solution with higher saturation in Fig. 2 (1st line of table I, the other solution would be indistinguishable on that 
same figure). It displays the result of the fit (full line) together with the contributions of the 1st (o^ 1 ), dotted line) 
and lst+2nd (a^ +a^ 2 \ dashed line) orders keeping the values of the parameters found for the full solution. We can 
see that these orders give distinguishable contributions. 



da 

dx 



dxdQ 2 dxjdki, 



1 We follow the same procedure as in Ref. |2£J, namely one HI point at k T > 5 GeV (7.3 1CT 4 ) and four ZEUS points (x = 4. 1CT 4 , 
and the three highest-x points) were not taken into account in the fits. The three highest— a; points for ZEUS cannot be described 
by a small— rr. approach probably because the x-value is too high (a; > 10 — 2 ). The other points cannot be de scrib ed because of large 
correction factors. Note that similar discrepancies appear also in other types of fitting procedures, e.g. in Ref. l25l . 



fit 


A 


Y 


MzEVS 


Mm 


X 2 (/dof) 


sat. 
weak sat. 


0.402 ± 0.036 ± 0.024 
0.370 ± 0.032 ± 0.015 


-0.82 ± 0.36 ± 0.01 
8.23 ± 0.48 ± 0.03 


34.3 ± 8.5 ± 7.0 
1136 ± 272 ± 2 


31.7 ± 8.4 ± 8.7 
1042 ± 238 ± 78 


6.8 (/ll) 
8.3 (/ll) 



TABLE I: Results of the fits to the Hl/ZEUS data for the GBW model. We find two independent solutions showing either 
significant (1st line) or only weak (2nd line) saturation parameters (see text). Statistical and systematic errors are indicated 
in that order for each parameter. 




FIG. 2: Results of the fit with significant saturation. Upper left: HI data kr > 3.5 GeV, upper right: HI data kr>5 GeV, lower 
left and right: ZEUS data in linear and logarithmic (showing the discrepancy at high- a;) scales. We display the result of the fit 
in full lines, together with the contribution of 1st and lst+2nd orders (resp. dotted and dashed lines). The cross-sections are 
given in nanobarns. 



Let us now discuss each solution. The first solution (1st line of Table I) shows a sizeably larger value of A than 
the i*2 parametrization ||| (cf. A^ 2 =0.288), while the value of Yq is found to be completely different (cf. Yq 2 =8.1), 
leading to more significant saturation effects. The second solution (2nd line of Table I) shows values of A and Yq 
more compatible with the F% result, even if A remains somewhat larger. The normalizations for HI and ZEUS data 
are found to be compatible in both cases. By contrast, the second and absolute x 2 minimum corresponds to fully 
significant (if not very big) saturation effects. To analyse these features more in detail, we study how the two x 2 
minima appear order by order in the expansion of the GBW cross-section using formula l|13Jl . In Table II, we show 



order 


A 


Y 


MzEVS 


Mux 


X 2 (/dof) 


1 

2 sat 
3 (sat.) 
3 (weak) 


0.372 ± 0.066 ± 0.015 
0.452 ± 0.027 ± 0.027 
0.447 ± 0.025 ± 0.026 
0.374 ± 0.028 ± 0.016 


-1.71 ± 0.47 ± 0.01 
-1.56 ± 0.33 ± 0.02 
6.85 ± 0.32 ± 0.12 


53.7 ± 10.9 ± 7.6 
20.2 ± 5.8 ± 4.9 
22.2 ± 3.0 ± 5.3 
693 ± 110 ± 78 


49.3 ± 11.5 ± 10.4' 
19.1 ± 5.2 ± 6.0 
20.9 ± 2.6 ± 6.6 

637 ± 94 ± 35 


8.3 (/12) 
6.8 (/ll) 
6.8 (/12) 
8.2 (/12) 



TABLE II: Results of the fits to the Hl/ZEUS data for the GBW models with different truncations. 1st line: 1st order (no 
saturation), 2nd line: 2nd order (~ GLR model), 3rd and 4th line: 3rd order with either significant or weak saturation effects. 



the results of the fits when we truncate the expansion up to first, second and third order. The first order 2 (1st line 
in Table II), which corresponds to the dilute limit ©, is very close to the 2nd minimum of the full model, showing 
that it indeed corresponds to weak saturation effects. This also confirms that a mere BFKL description [2£j of da/dx, 
similar to the first order approximation ©, fits well the same data. At second order, which corresponds to the GLR 
version of the model, only one minimum appears and it is closer to the first minimum of Table I. At third order, one 
recovers the two minima observed with the full model. It is clear that saturation corrections are needed to see the first 
minima appear and this confirms that it involves significant saturation effects. The weak saturation solution should 
a priori be present regardless of the number of terms in the expansion, yet it is not the case when the truncation is 
done at second order. Our interpretation of this feature is that the cross-section (0} is sensitive to large dipoles and 
to the high— r e ff behavior of the truncated dipole-dipole cross-section even when the saturation radius Rq is large; 
this is due to the fluctuations of the jet dipole distribution (JHJ around its mean value. 

Let us make some additionnal comments, (i) We used ot 2 s (in the overall factor) running at one loop, and for 
completion, we checked the results of the fits while keeping as constant. We find the same quality of the fits but 
with the effective saturation intercept enhanced by about 30% for the first solution, (ii) We have also studied the 
other models for effective radii proposed in which lead to fits of the same quality with comparable values for 
the paramaters. (iii) The physical normalization (see equations (5) and (9)) is A r exp(— XYq). We checked that we 
obtain a set of consistent values around 50 for this physical normalisation for all fits. The disperse values of TV are 
compensated by the values of Y . 

In Fig 3. we plot the different saturation scales as a function of the physical rapidity interval Y. We also display the 
saturation scale parametrization obtained by the GBW fit of the proton structure function F 2 |8| . The weak solution 
is compatible with the one from F% in the physical range from five to ten units in rapidity, the rapidity dependance 
is however somewhat stronger. This solution favors a universal saturation scale regardless of the presence of a soft 
scale in the problem, as proposed in |12J where 7*7* data are described with the saturation scale obtained from F 2 . 
The strong saturation solution gives a curve that lies significantly above the weak and F2 curves, indicating that in 
forward jets, the saturation scale is different and higher than for deep inelastic scattering, as claimed in |ll|. 

At this point it is important to compare the range of Q s = 1/i? 2 , with the experimental transverse momentum cuts. 
Notice that these cuts (3.5 and 5 GeV) lie approximately in between the two possible saturation scales. It is thus 
clear why the saturation effects are different for both solutions. 



IV. PREDICTIONS FOR THE TEVATRON AND THE LHC 



The fact that the saturation corrections can be sizeable for two-scale processes suggests that the predictions for 
related hard cross-sections at hadronic colliders will show striking differences with the weak saturation case. Let us for 
instance consider the Mueller-Navelet jet production at Tevatron and LHC following the approach of 0]. Indeed, 
the Mueller-Navelet jets are directly related to our discussion and, at the level of the two-scale hard cross-section 
amounts to replacing the photon probe by another forward jet. One writes 



o-jj(h 



T1,KT2 



Ary) 



16 



d 2 r d 2 f <f> J (r,]4i) <t> J (f,k^ 2 ) 



1 — exp 



r lff( r > r ) 
4i?2(Ary) 



(15) 



2 In this case, there is one less parameter as e appears only in the normalization and is absorbed in A/". 



2 4 6 8 10 12 

Y 



FIG. 3: Saturation scales log Q%/Qo = A(Y — Yq). We display the saturation scales for both solutions of our fits (full lines with 
errors indicated by dashed lines) and compare with the result for F2 (dotted line). The weak saturation scale and the one for 
F2 are compatible: lower curves. The other solution lies above, yielding stronger saturation effects. 



where Ary is the rapidity interval between the two jets, and kxi and fc<T2 are their lowest transverse momenta. We 
thus consider cross-sections integrated over the transverse momenta of the jets with lower bounds kn and kx2- In 
order to appreciate more quantitatively the influence of saturation, it is more convenient to consider the quantities 
Tli/j defined as 

_ <j{k T i,k T 2, Arjj) 
a{k T1 ,k T2 ,AT] j ) 

i.e. the cross-section ratios for two different values of the rapidity interval. These ratios display in a clear way 
the saturation effects, they also correspond to possible experimental observables if one changes the center-of-mass 
energy since they can be obtained from measurements at fixed values of the jet light-cone momentum and thus are 
independent of the parton densities the incident hadrons. Measuring two values of the cross-section for identical jet 
kinematics and different rapidity intervals A 7^ and Arjj and dividing them gives the ratio TZ-i/j- Such observables 
have actually been used for a study of Mueller-Navelet jets for testing BFKL predictions at the Tevatron pfj l2fil|. 
However, it is known that hadronization effects play a role in this kind of measurements and have to be taken into 
account |27| . 

In Fig 4, we show the resulting ratios, when kri = kx2 = kx, for 72-5/2 which corresponds to accessible rapidity 
intervals at the Tevatron and TZg/2 which corresponds to realistic rapidity intervals for the LHC. The curves are for 
both saturation solutions together with the prediction obtained from the GBW parametrization of F2. At high scales, 
both models lead to similar values of 7Z, larger than the ones extracted 0] from F 2 . This reflects the higher saturation 
intercept systematically found with forward jets in both cases. 

A clear difference between the two options (saturation being weak or not) appear in the drop due to saturation 
which occurs for kx^Qs- For both Tevatron and LHC kinematical configurations, the weak saturation solution would 
only give effects at rather and probably too low values of fc-r to be observable. By contrast, the stronger saturation 
solution would give rise to saturation effects for kr of order 3 GeV for the Tevatron and 30 GeV for the LHC. While 
the effects could be marginal at Tevatron, the larger value of fe T for the LHC provides a way to distinguish between 
the two solutions and to test the models. Hence the observation or not of saturation effects for Mueller-Navelet jets 
at LHC is a promising discriminant between the two options. 
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FIG. 4: Cross-section ratios TZi/j. R5/2 and Rs/2 are defined for rapidity intervals 2 and 5 for the Tevatron (3 lower curves), 
and 2 and 8 for the LHC (3 upper curves). The full lines correspond to the significant saturation solution, first line of Table 
I. The dashed lines are for the weak saturation solution, second line of Table I. The dotted lines correspond to the saturation 
parametrization of F2. 



V. CONCLUSION 



Let us briefly summarize the main results of our study. We described the published HI and ZEUS forward-jet data 
using a saturation model based on the GBW formalism. We find two possible fits to the data with either significant 
or weak saturation. The weak saturation solution would be in favor of a universal QCD saturation scale. The solution 
with significant saturation effects would prove process dependent saturation scales. 

Hence the confirmation or not of two different saturation scales, one for one-hard-scale processes (e.g. for F 2 ), and 
one for two-hard-scale processes (e.g. for forward-jets) would suggest interesting theoretical questions, for instance 
whether the saturation scale is universaly related e.g. to Aqcd or could be related to the hard probes initiating the 
process. 

Both fits lead to different predictions at the Tevatron, and even more different at the LHC where the effect is 
sizeable for larger jet transverse momentum. The measurement of Ti-i/j would imply running the accelerators at 
different center-of-mass energies. It would allow to test if saturation is stronger for harder processes compared to deep 
inelastic scattering. It will also be interesting to test the two solutions using the recently presented preliminary HI 
and ZEUS forward-jet and forward-pion data, when available. 
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